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Abstract 

We present exact calculations of the Potts model partition function Z(G,q,v) for ar- 
bitrary q and temperature-like variable v on n-vertex square-lattice strip graphs G for a 
variety of transverse widths L t and for arbitrarily great length Lg, with free longitudinal 
boundary conditions and free and periodic transverse boundary conditions. These have the 
form Z(G,q,v) = Y^=i' X c z,G,j{^z,G,j) Ll - We give general formulas for Nz,gj and its spe- 
cialization to v = — 1 for arbitrary L t for both types of boundary conditions, as well as other 
general structural results on Z. The free energy is calculated exactly for the infinite-length 
limit of the graphs, and the thermodynamics is discussed. It is shown how the internal energy 
calculated for the case of cylindrical boundary conditions is connected with critical quanti- 
ties for the Potts model on the infinite square lattice. Considering the full generalization to 
arbitrary complex q and v, we determine the singular locus £>, arising as the accumulation 
set of partition function zeros as Lg — > oo, in the q plane for fixed v and in the v plane for 
fixed q. 
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1 Introduction 



The g-state Potts model has served as a valuable model for the study of phase transitions 
and critical phenomena [0-0. On a lattice, or, more generally, on a (connected) graph G, 
at temperature T, this model is defined by the partition function 

Z(G,q,v) = J2^ m (1-1) 

{cr„} 

with the (zero-field) Hamiltonian 

n = -jJ2^ (i.2) 

where <Ji = 1, q are the spin variables on each vertex % 6 G; (3 — (A^T) -1 ; and (ij) denotes 
pairs of adjacent vertices. The graph G = G(V, E) is defined by its vertex set V and its edge 
set E\ we denote the number of vertices of G as n = n(G) = \V\ and the number of edges of 
G as e(G) — \E\. We use the notation 

K = (3 J , a = e K = u- 1 , v = a-l (1.3) 

so that the physical ranges are (i) v > corresponding to oo > T > for the Potts ferro- 
magnet, and (ii) — 1 < v < 0, corresponding to < T < oo for the Potts antiferromagnet. 
One defines the (reduced) free energy per site / = —j3F, where F is the actual free energy, 
via 

f(G,q,v) = n- 1 \n[Z(G,q,v)] (1.4) 

and, in the n — » oo limit, 

f({G},q,v) = lim n- 1 ln[Z(G,g,^)] (1.5) 

where we use the symbol {G} to denote lim n ^ 00 G for a given family of graphs. 

Let G' = (V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex 
set V and an edge set E' C E 1 . Then Z(G, g, f ) can be written as the sum || 

Z(G,q,v)= q k{G ' ] v e{G,) (1.6) 

G'CG 

where k(G') and e(G') denote respectively the number of connected components and the 
number of edges of G'. Since we only consider connected graphs G, we have k(G) = 1. The 
formula (|1.6|) enables one to generalize q from Z + to M + and, indeed, to C. The Potts model 
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partition function on a graph G is essentially equivalent to the Tutte polynomial ||-f7j and 
Whitney rank polynomial @, f8fl-||10|| for this graph. 



In this paper we shall present exact calculations of the Potts model partition function 
Z(G,q,v), for arbitrary q and v, on square-lattice strip graphs G of various widths L t and 
arbitrarily great length Lg, with boundary conditions that are (i) free in both longitudinal 
and transverse directions, denoted "open" and (ii) periodic in the transverse direction and 
free in the longitudinal direction, denoted "cylindrical" . Specifically, we shall consider strips 
with L t vertices along a transverse slice (hence L t edges on this slice in the case of cylindrical 
boundary conditions and L t — 1 edges in the case of free boundary conditions) and L e vertices 
along a longitudinal slice. We shall often denote the strip with open boundary conditions as 
(L t )p x (Li)-p and the strip with cylindrical boundary conditions as (L t )p x (L^p. Previous 
related calculations of Z(G, q, v) for arbitrary q and v on fixed-width, arbitrary-length lattice 
strip graphs are in |ll]]-|[T3| for the square lattice, in |14| for the triangular lattice, in (T^j for 



the honeycomb lattice, and in [16| for the square lattice with next-nearest-neighbor edges 
(bonds). Some general properties are given in [17, 18]. Calculations of Z(G, q, v) for arbitrary 
q and v on finite sections of lattices have been given in []19 , p0| and used for studies of zeros 
in the q and v plane (the latter being complex-temperature or Fisher zeros 0]); for fixed 



positive integral values of q, such calculations have been given, e.g., in [p2]|-[p8~|. Since the 
Ising case q = 2 is exactly solvable on two-dimensional lattices, one can calculate exactly 
the complex-temperature phase boundaries that are the accumulation sets of the complex- 
temperature zeros in the thermodynamic limit (e.g., [29]-[41] and references therein). 

There are several motivations for the present work. Clearly, new exact calculations of 
Potts model partition functions are of value in their own right. Since the results apply for 
arbitrary length, one can take the limit of infinite length and calculate the free energy and 
specific heat for the Potts model with any q and various widths. In addition, it was shown 
|12|1 that these calculations can give insight into the complex-temperature phase diagram of 
the 2D Potts model on the given lattice. This is useful since the 2D Potts model has never 
been solved except in the (zero-field) q = 2 Ising case. 

Various special cases of the Potts model partition function are of interest. One special 
case is the zero-temperature limit of the Potts antiferromagnet (AF), i.e., v = — 1. For 
sufficiently large q, on a given lattice or graph G, this exhibits nonzero ground state entropy 
(without frustration). Such nonzero ground state entropy is an exception to the third law 
of thermodynamics f42| , 43]. It is equivalent to a ground state degeneracy per site (vertex), 
W > 1, since So = fcelniy. The T = (i.e., v = —1) partition function of the above- 
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mentioned g-state Potts antiferromagnet (PAF) on G satisfies 



Z(G,q,-l)=P(G,q) (1.7) 

where P(G, q) is the chromatic polynomial (in q) expressing the number of ways of coloring 
the vertices of the graph G with q colors such that no two adjacent vertices have the same 



color [0, 44, 45 . The minimum number of colors necessary for this coloring is the chromatic 



number of G, denoted x(G). We have 

W{{G},q) = lim P(G,q) 1/n . (1.8 

n— >oo 

Some previous related works on chromatic polynomials and the W function are [46] -[f77| 



Using the formula ( |1.6|) for Z(G, q, v), one can generalize q from Z + not just to M + but to C 
and a from its physical ferromagnetic and antiferromagnetic ranges < v < oo and — 1 < 
v < to v G C. A subset of the zeros of Z in the two-complex dimensional space C 2 defined 
by the pair of variables (q, v) can form an accumulation set in the n — > oo limit, denoted B, 
which is the continuous locus of points where the free energy is nonanalytic. This locus is 
determined as the solution to a certain {G}-dependent equation fll], [OH- For a given value 
of v, one can consider this locus in the q plane, and we denote it as B q ({G}, v). In the special 
case v — — 1 where the partition function is equal to the chromatic polynomial, the zeros in 
q are the chromatic zeros, and B q ({G},v = —1) is their continuous accumulation set in the 
n — > oo limit. Previous works have given exact calculations of the chromatic polynomials and 
nonanalytic loci B q for various families of graphs. With the exact Potts partition function, 
one can study B q for arbitrary temperature in both the antiferromagnetic and ferromagnetic 
cases and, for a given value of q, one can study the continuous accumulation set of the zeros 
of Z(G, q, v) in the v plane; this will be denoted B V ({G}, q). (Where the context is clear, we 
shall refer to these accumulation sets in the q and v plane as just B.) 

Following the notation in and our other earlier works on B q ({G}) for v = —1, we 
denote the maximal region in the complex q plane to which one can analytically continue 
the function ^({G}, q) from physical values where there is nonzero ground state entropy as 
R\. The maximal value of q where B q intersects the (positive) real axis was labelled q c ({G}). 
Thus, region Ri includes the positive real axis for q > q c ({G}). Correspondingly, in the stud- 
ies of complex-temperature phase diagrams of spin models a nomenclature was used 
according to which the complex-temperature extension(CTE) of the physical paramagnetic 
phase was denoted as (CTE)PM, which will simply be labelled PM here, the extension being 
understood, and similarly with ferromagnetic (FM) and antiferromagnetic (AFM) phases; 
other complex-temperature phases, having no overlap with any physical phase, were denoted 
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Oj (for "other"), with j indexing the particular phase. Here we shall continue to use this 
notation for the respective slices of B in the q and v planes. 

We record some special values of Z(G, q, v) below, beginning with the q = special case 

Z(G,0,v) = (1.9) 

which implies that Z(G, q, v) has an overall factor of q, as is also clear from ( |1.6|) . In general 
(and for all the graphs considered here), this is the only overall factor that it has. We also 
have 

Z(G, l,v)= V<G>) = a<G) ■ ( L1 °) 

G'CG 

For temperature T = oo, i.e., v = 0, 

Z(G, q, 0) = <f (G) . (1.11) 
where n(G) is the total number of sites of G. Further, 

X(G)-1 

Z(G, q, -1) = P(G, ?) = [ J] («- s )] ^( G > ^ (1-12) 

where g) is a polynomial in g of degree n(G) — 

Another basic property, evident from eq. (|1.6| ), is that (i) the zeros of Z(G, q, v) in q for 
real v and hence also the continuous accumulation set B q are invariant under the complex 
conjugation q —>■ q*; (ii) the zeros of Z(G, q, v) in v or equivalently a for real q and hence also 
the continuous accumulation set B v are invariant under the complex conjugation v — > v*. In 
certain cases one must take account of the noncommutativity ]5(| 

lim lim Z{G,q,v) 1/n ^ lim lim Z(G,g,w) 1/n . (1.13) 

n— >oo 5— >q s q-^Qs n— >oo 

A general form for the Potts model partition function for the strip graphs considered 
here, or more generally, for recursively defined families of graphs comprised of m repeated 
subunits (e.g. the columns of squares of transverse width L t vertices that are repeated Le 
times in the longitudinal direction to form an Lt x L t strip of a regular lattice with some 
specified boundary conditions), is [ |T^ 

Z(G,q,v)= <*A*Gj) m (1-14) 

i=i 

where the coefficients cqj and corresponding terms Xgj, as well as the total number Nz,g,x 
of these terms, depend on the type of recursive graph G (width and boundary conditions) 
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but not its length. In the special case v = — 1 where Z reduces to the chromatic polynomial 
(zero-temperature Potts antiferro magnet), eq. ( 1.14Q reduces to the form p2| 



P(G,q)= cgj{\pmT. (1.15) 

3=1 

For the lattice strips of interest here, we define the following explicit notation. Let N z , S q,BCt BC t> L t ,x 
denote the total number of A's for the square-lattice strip with the transverse (t) and longitu- 
dinal (£) boundary conditions BC t and BCi of width L t . Henceforth where no confusion will 
result, we shall suppress the A subscript. The explicit labels are Nz, S q,FF,L t an d Nz,tH,FF,L t f° r 
the strips of the square and triangular lattices with free boundary conditions, and N z , S q,PF,L t 
and Nz,tri,PF,L t f° r the strips of these respective lattices with cylindrical boundary condi- 
tions. In the special case v = — 1, our nomenclature for the corresponding total numbers 
of A's in the chromatic polynomials are listed below, together with the labels used for these 



numbers in 17 



N P , sq ,FF,L t = SqFree(L t ) (1.16) 

N P . trimt = TriFree(Li) (1.17) 

N P ^ PFM = SqCyl(L t ) (1.18) 

N PM:PFM ee TriCyl(L t ) . (1.19) 

The Potts ferromagnet has a zero-temperature phase transition in the Li — >• oo limit of 
the strip graphs considered here. However, in contrast to the infinite-length, finite-width 
lattice strips with periodic (or twisted periodic) longitudinal boundary conditions, where 
this property was reflected in the feature that B in the u = 1/a plane passes through the 
T = point u = and hence is noncompact in the a plane, for infinite-length strips with 
free longitudinal boundary conditions, B is generically compact in the a (equivalently, the 
v) plane. 

2 General Structural Properties 

In this section we give some results for the numbers of A's that enter in the respective 
chromatic polynomials and Potts model partition functions for the strips of interest here. 



Recall the definitions of the Catalan number C n and Motzkin number M n ||78||-f8T 



c«=4t( 2 ;) (2.D 

n + 1 \ n 
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M„ 



i=o KJ J 



(2.2) 



These numbers occur in many combinatoric applications ||78|]-]B3[]. Among these is the con- 
struction of non-intersecting chords on a circle; the number of ways of connecting a subset of 
n points on a circle by non-intersecting chords is M n , while the number of ways of completely 
connecting 2n points on the circle by non-intersecting chords is C n . Summing over subsets 
of points connected by chords, this yields the known relation 



E 

A:=0 



n 
2k 



(2.3) 



where [u] denotes the integral part of v. These numbers have the asymptotic behaviors 



l + 0(n' 1 ) 



as n 



oo 



M„ 



3 3/2 2 -l 7r -l/2 n -3/2 3 n 



1 + 0(n 



as n 



oo 



(2.4) 



(2.5) 



Ref. |18| derived expressions for the total number of A's in Z and P with periodic bound- 
ary conditions in the longitudinal direction and free boundary conditions in the transverse 
direction. We denote such quantities respectively as N ZtSq tri,FP,L t an d Np jSqtrijFPtLt , where 
sqtri means either square or triangular: 



z ,sqtri,FP,Lt 



2L t 



(2.6) 



l 2 



Np :S qtri,FP,Lt ~ ^{Lt ~ l)! 



(Lt-j) 



3=0 



(j\y(L t -2j)\ ■ 



These have the leading asymptotic behaviors Nz )Sq tri,FP,L t ~ L t 4 Lt and Np tSqtTi ^FP,L 
L t 3 Lt as L t —y oo. 



(2.7) 



2.1 Np iGiB c,L t 

Recall that for the open strip of the triangular lattice of width L t , Ref. [17] obtained the 
result 

N P ,tri,FF,L t = M Lt _! . (2.1.1) 

The proof of this was based on the property that there is a one-one correspondence between 
the number of distinct eigenvalues A of the transfer matrix and the number of non-crossing, 
non- nearest-neighbor partitions of the set {1, 2, L t }. But the latter number is precisely 
the Motzkin number Ml \. 
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Table 1: Numbers of A's for the chromatic polynomial for the strips of the triangular and square lattices 
having free and cylindrical boundary conditions and various widths L t . 





Np,tri,FF,L t 


Np, S q,FF,L t 


— ^P,tri,FF,L t 


d Lt 


N^P,tri,PF,L t 


Np iSq ,pF,L t 


-Np it ri,PF,L t 


Lt^P,tri,PF,L t 

-d Lt 


1 


1 


1 


1 


1 


1 


1 


1 





2 


1 


1 


1 


L 1 


1 


1 


1 


1 


3 


2 


2 


2 


1 


1 


1 


1 


2 


4 


4 


3 


2 


3 


2 


2 


2 


5 


5 


9 


7 


5 


6 


2 


2 


2 


4 


6 


21 


13 


5 


15 


5 


5 


5 


15 


7 


51 


32 


13 


36 


6 


6 


6 


6 


8 


127 


70 


13 


91 


15 


14 


13 


29 


9 


323 


179 


35 


232 


28 


22 


16 


20 


10 


835 


435 


35 


603 


67 


51 


35 


67 


11 


2188 


1142 


96 


1585 


145 


95 


45 


10 


12 


5798 


2947 


96 


4213 


368 


232 


96 


203 


13 


15511 


7889 


267 


11298 


870 


498 


126 


12 


14 


41835 


21051 


267 


30537 


2211 


1239 


267 


417 



Now the number of A's for the open square-lattice strip with a given width L t is less than 
for the open triangular strip of this width because some of the symmetry under reflection 
about the longitudinal axis, which renders some of the partitions equivalent to each other. We 
proceed to derive this number. Because of this reflection symmetry 2Np^ FFyFt ~ Npjri,FF,L t 
gives the number of non-crossing non-nearest-neighbor partitions for a transverse slice that 
is symmetric under this reflection. 

Using this observation, we proceed to derive the following result: 

Theorem 1 ^ 

2N PiSq>FF>Lt - N P:triiFF:Lt = -N PsqtH Fn L i + 1] . (2.1.2) 

Proof. Let us denote 2N PjSg:FF:Lt — N PjtriiFFjLt by X Lt for short. To show the general 
method of the proof, we first list, for small L t , the sets Px Lt of partitions having the above- 
mentioned reflection symmetry: Px 1 = {1}, Px 2 — {!}, Px 3 = {1,^1,3}, Px 4 = {l,#i,4}, 

Px 5 = {M2,4 A )5 , 5 2 ,4<Sl,5 A, 3 ,5}, P X 6 = {1, $2,5, #1,6, S 2 , 5 Sl,6, <5l,3&M>}- Here ll Can be Seen 

that going from the partition for an odd width L t to the partition for the next greater width, 
L t + 1, one only has to split the central vertex into two vertices (one connected only to 
vertices left of center, the other only to vertices right of center) without otherwise changing 
the partitions. Hence, 

X Lt+1 =X Lt for odd L t . (2.1.3) 

Going from the partition for an even width L t to the partition for the next greater L t + 1, 
there are two possibilities. The first possibility is that there is no delta function connecting 
any vertex in the set {1, 2, .., L t /2} and the corresponding vertex in the set {L t /2 + 1, L t /2 + 
2, ..,L t }. Let us denote the number of these kinds of partitions by X LtA , which is equal to 
N R 

tri,FF,L t /2- in this case, one can add a vertex in the middle. The second possibility is that 
there is more than one delta function connecting the two vertex sets mentioned above. Let 
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us denote the number of these kinds of partitions by Xi, t f,- Clearly, X^ ua + Xi, u b = Xp t for 
the case of even L t considered here. In this case, one can either add a vertex in the middle 
without any further delta function or have one more delta function connecting this middle 
vertex and the nearest vertex pair which are already connected by a delta function. 

Another way to get a partition for L t + 1 vertices is to start with a basis consisting of 
L t — 1 vertices. Here one splits the central vertex for L t — 1 into three vertices with the upper 
and lower vertices connected by a delta function. The number of these types of partition is 
X Lt _ x . With eqs. ( pXlp and ( gXp , we obtain 



X 



Lt+l 



Xh t ,a + ^X Luh + X Lt _i 

M±_ x + 2{X Lt _ 1 - M il _ 1 ) + X Lt 



2 1 



for even L 



t ■ 



Setting L t = 2£ in eq. (|2.1.6| ), we have 

X 2e+1 = 3 £ - ??- l M - 3^ 2 Mx 



21+2 



(2.1.4) 
(2.1.5) 
(2.1.6) 



(2.1.7) 



We observe from eq. ( |2.1.6 ) that 2X 2 e+i has the same recursion relation and initial values 
as iVpL t+ljA in eq. (5.1) in [Pfj , which proves eq. ( |2.1.2| ). □ 



We remark that this proof evidently links together results for two different sets of longitudinal 
boundary conditions. 



Theorem 2 



hr, (H-j) 



N PjSq>FFjLt =-M Lt „ l + -(L r t -l)\J2 



3=0 



where 



Lt + l 



[2.1.8) 



(2.1.9) 



Proof This follows from Theorem 1 and the expressions for Npj ri ^F,L t an d Np jSqtri)F p^ Lt+1 y 2 ]- 



□ . 



Corollary 1. N PtSq)FFjLt has the asymptotic behavior 



N. 



P,sq,FF,L t 



31/23-2^-1/2^-3/23^ h + o( L -i 



as 



(2.1.10) 



Proof. This follows directly from the asymptotic behaviors given above for the Motzkin 
number and for N PtSqtT ^ FP ^ Lt) which show that as L t — > 00, the third term, N PjSqtritFP ^ Lt+ iy 2 ], 



is negligibly small compared with the first two in eq. (|2.1.2| ). The latter, in turn, is a 
consequence of the fact that the leading exponential asymptotic behavior, as L t — > 00, of 
the dominant terms is ~ 3 Lt while that of the third term is ~ 3 it//2 . 
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Corollary 2. 

lim = I . (2.1.11) 

Nt-nxi Npjri,FF,N t 2 

Proof. This follows from ( |2.1.2| ) and the asymptotic behaviors given above. 

In |17J it was conjectured that lmij^oo SqFree(L 4 )/TriFree(L t ) = 1/2. Our theorem and 
corollary 2 prove this conjecture. 

We next consider the strips with cylindrical boundary conditions and the quantity 2Np tSqt PF,L t 
N P . tri,PF,L t - We find that the following relation holds for 1 < L t < 14 and conjecture that it 
holds for all widths L t : 

Conjecture 1 For arbitrary L t , 

\ \ N p sq triFP± for even ^ 

2N PM - *W«. = < iNp ^ p ' - iRh p for odd U > 3 (2 - L12) 



where R n is the Riordan number, which may be defined via the generating function ||81||,||84 

v ' n=0 

or, for n > 2, by |35| 

n-l 

R n = (-1)™ ^{-ly^Mj for n > 2 . (2.1.14) 

i=i 

Note that Ri t is equal to the quantity di t given in Table [l] for L t > 2. 

Theorem 3. d L 1 

Np t tH,PF,L t = — — for prime L t . (2.1.15) 



Proof As shown in JT7| , d^ t is the number of non-crossing non- nearest-neighbor partitions 
of the vertex set {1, ...,L t } with periodic boundary conditions on the set. Now, when L t is 
a prime number, all partitions except the partition 1 (i.e., all blocks being singletons) have 
an L r fold degeneracy. Therefore, 

LtNp,tri,PF,L t - d Lt = L t - 1 for prime L t (2.1.16) 

and the theorem follows. □ 

In |17| it was conjectured that the number TriCyl(L 4 ), i.e., Np }tr i,PF,L t , behaves asymp- 
totically like di t /L t as L t — > oo. Our Theorem 3 ( |2.1.15| ) proves this conjecture for prime 
L t . 

From eqs. ( 2.1.1 2| ) and Q2.1.15Q , we infer 
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Table 2: Numbers of A's for the Potts model partition function for the strips of the triangular and square 
lattices having free and cylindrical boundary conditions and various widths L t . 



L t 


Nz,tri,FF,L t 


Nz,sq,FF,L t 


2Nz,sq,FF,L t 
—Nz,tri,FF,L t 


Nz,tri,PF,L t 


Nz,sq,PF,L t 


2N Z ,sq,PF,L t 
—Nz,tri,PF,L t 


LtNz,tri,PF,L t 
—Nz,tri,FF,L t 


1 


1 


1 


1 


1 


1 


1 





2 


2 


2 


2 


2 


2 


2 


2 


3 


5 


4 


3 


3 


3 


3 


4 


4 


14 


10 


6 


6 


6 


6 


10 


5 


42 


26 


10 


10 


10 


10 


8 


6 


132 


76 


20 


28 


24 


20 


36 


7 


429 


232 


35 


63 


49 


35 


12 


8 


1430 


750 


70 


190 


130 


70 


90 


9 


4862 


2494 


126 


546 


336 


126 


52 


10 


16796 


8524 


252 


1708 


980 


252 


284 


11 


58786 


29624 


462 


5346 


2904 


462 


20 


12 


208012 


104468 


924 


17428 


9176 


924 


1124 


13 


742900 


372308 


1716 


57148 


29432 


1716 


24 



Conjecture 2 For prime L t > 3, 

Np,q,PF, Lt = \ ( $h±£z± + \N psqtn Fp ^ - l -R^) . (2.1.17) 



Note that Conjecture 2 is consistent with the conjecture made in |I7| that in the limit as 
L t — > oo the ratio SqCyl(L t )/TriCyl(L t ) = Np tS g t PF t L t /Np ttr i,PF,L t is equal to 1/2. For 
non-prime L t , the right-hand sides of eqs. ( |2.1.15|) and ( |2.1.17|) are the lower bounds for 



N P . tri,PF,L t and Np tS q : PF : L t , respectively. However, we have not yet obtained formulas for 
N Pt tri,PF,L t or N P)S g jPF>Lt for arbitrary L t . 

2.2 N Z ,G,BC,L t 

Here we carry out a similar analysis for the total number of A's in the Potts model partition 
function or equivalently, the Tutte polynomial for various lattice strips. First, we list as the 
left-most column in Table [| the number of A's that enter in this partition function for the 
open strip of the triangular lattice of width L t , i.e., Nz,tH,FF,L t - This number is equal to 
the number of non-crossing partitions of the set {1, 2, L t }, which is precisely the Catalan 
number. From this observations, it follows that |T^J 

N Z ,tri,FF,L t = C Lt . (2.2.1) 

In the zero-temperature limit for the antiferromagnetic case, one can work in a subspace that 



corresponds, after a suitable change of basis, to only non-nearest-neighbor partitions [17 . 
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In this limit, the number of A's is thus reduced to the non-crossing, non-nearest-neighbor 
partitions given by M Lt _i in Q2.1.1J) . 

Just as was the case for the chromatic polynomial, the number of A's for the open square- 
lattice strip with a given width L t is less than for the open triangular strip of this width 
because of the symmetry under reflection about the longitudinal axis, which renders some 
of the partitions equivalent to each other. We proceed to derive this number. Because 
2N z , S q,FF,Lt ~ Nz,tri,FF,L t gives the number of non-crossing partitions for a slice of the trans- 
verse vertices which is symmetric under the reflection symmetry, we have 

Theorem 4 

/ N Z sqtri FPM eVeI1 Lt 

2N ZtSq , FF>Lt - N WFM = t i ' ' ' 2 (2.2.2) 

I 2 Z,sqtri,FP, ' 1 



where N z , S qtri,FP,L t was listed in eq. ( j2.6|) from |T8 |. 

Proof Let us denote 2Nz, S q,FF,L t — Nz,tri,FF,L t by Yi t for simplicity, and list, for small 
L t , the sets Py i( of partitions having reflection symmetry: Py x = {1}, Py 2 = {1, #1,2}, Py 3 = 

{1,^1,3,^1,2,3}, Py 4 — {l,^2,3,^l,4,^2,3^1,4,5l,253,4,^l,2,3,4,5}, Py 5 = {1,^2,4,^1,5,^2,4^1,5,^1,2^4,5, 

^2,3,4, ^1,2,4,5, ^2,3,4^1,5, £1,3,5, £1, 2,3,4,5}- Here it can be seen that going from an odd L t to L t + 1, 
the number of partitions is doubled because one can split the central vertex to two vertices 
with or without a delta function for these two vertices. Therefore, 

Y Lt+1 = 2Y Lt for odd L t . (2.2.3) 

Going from the partition for an even L t to the partition for L t + 1, there are two possibil- 
ities. The first occurs if there is no delta function connecting any vertex in the vertex set 
{1,2,.., L t /2} and the corresponding vertex in vertex set {L t /2 + 1, L t /2 + 2, .., L t }. Let us 
denote the number of these kinds of partitions by Yi t:a , which is equal to Nz,tri,FF,L t /2- 

In 

this case, one can add a vertex in the middle. The second possibility occurs if there is more 
than one delta function connecting the two vertex sets mentioned above. Let us denote the 
number of this kind of partitions by Y Lt>b . Clearly Y Lua + Y Lub = Y Lt for the case of even L t 
considered here. For this case one can either add a vertex in the middle without any further 
delta function or have one more delta function connecting this middle vertex and the nearest 
vertex pair which are already connected by a delta function. Using eqs. ( |2.2.1| ) and ( |2.2.3| ), 
we obtain 

Y Lt+1 = Y Lua + 2Y Lub (2.2.4) 
= CL t +2{2Y Lt - 1 -CL t ) (2.2.5) 

2 2 

= 4F Lt _i - Cl± for even L t . (2.2.6) 

2 

Set L t = 2£ in eq. ( |2.2.6| ), we have 

y 2m = 4 l -4^ 1 C 1 -4^ 2 C 2 -...-Q (2.2.7) 

1 (21 + 2 \ . . 

2 1 J 

(2.2.9) 
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and 



21: 



(2.2.10) 



Using (|2.q) ) this completes the theorem. □ 

For even L t , N z sqtri Fp l± is equal to the entry in the central column of Pascal's triangle 
(with rows r > and elements on the r'th row given by rj for < s < r). For odd L t , 



(1/2)N 



Z,sqtri,FP, 



L f +i is equal to the entry in the next-to-central column of this triangle [| 



Theorem 5 



N 



Z,sq,FF,L t 



c Lt + ( L 4) 

2 



c, 



1 (Lt+1\ 



for even L t 
for odd L t 



(2.2.11) 



Proof The theorem follows from Theorem 4 and eq. ( 2 . 2 . 1| ) . □. 



The sequence of numbers Nz tS q,FF,L t is equal to the number N bT)Lt+ i of distinct bracelets 
(i.e. distinct modulo rotations and reflections) with L t + 1 black beads and L t white beads 
(given as sequence A007123 in ]S3| and previously as sequence M1218 in |S2fl). This can 
be understood as follows. The number of bracelets enumerated without dividing out by 
equivalence classes of reflections, i.e. the number that are distinct up to rotations only, with 
L t + 1 black beads and L t white beads, is given by the Catalan number C^- Of these, 
the number that are symmetric under reflection is (lJ) for even L t and ( l^U ) for odd L t . 

2 2 

Therefore, the number of distinct bracelets is the same as the sequence of numbers Nz, sq ,FF,L t - 
Corollary 3 

The number Nz, S q,FF,L t has the leading asymptotic behavior 



N 



Z,sq,FF,L t 



-n- l l 2 Lf' 2 k L * + 2 l ' 2 L t 2 L t + 0(- 



Proof This follows from theorem 5 and eq. 
Corollary 4 



hm N J^ FF ^ 

L t ->oo Nz,tri,FF,L t 



V 

!4j). Hence, 



as L t 



(2.2.12) 



(2.2.13) 



Proof This follows from ( |2.2.1| ) and our new result ( |2.2.11| ). 

Similarly, consider the quantity 2Nz, S q,PF,L t — Nz,tri,PF,L t f° r the strips with cylindrical 
boundary conditions. We find that this is the same as 2N ZjS q,FF,L t ~ Nz,tri,FF,L t f° r 1 ^ L t < 
13 and infer the following conjecture: 

Conjecture 3 For arbitrary L t , 



2N Z ,sq,PF,L t -N z ,tri,PF,L t 



Nz,sqtri,FP,±± 

In 

2 iV Z,sqtri,FP,±±±± 



for even L t 
for odd L t 



(2.2.14) 
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Theorem 6 



N 



Z,tri,PF,L t 



C Lt + 2(L t -l) 



for prime L 



t ■ 



(2.2.15) 



Proof As shown in | l7| , Nz t tH,FF,L t is the number of non-crossing partitions of vertex set 
{!,... ,L t }. Now if L t is a prime number, except for the partitions 1 (i.e., all blocks being 
singletons) and 8\ t 2,...,L t (i- e -> a unique block), there is an L r fold degeneracy for the rest of 
the partitions. Therefore, 



LtN Z,tri,PF,L t 

and the theorem follows. □ 



N 



Z,tri,FF,L t 



2(L t 



for prime L t 



(2.2.16) 



From reference to f83| , we observe that the numbers N Z jri,PF,L t that we have calculated 
are equal to the sequence listed as A054357 in [53] which enumerates the number of tree 
graphs with L t edges modulo equivalence under 2-coloring (so that, for example, this enu- 
meration includes only one path graph if L t is odd and includes two path graphs if L t is 
even). We also observe that our result ( [2.2.15|) is the special case, for prime L t , of this 
sequence. A general analytic formula for the above sequence is given in [|55 ]. Based on these 
connections, we conjecture that 



Conjecture 4 For arbitrary L t 

N z 

,tri,PF,Lt 



1 



d\L t ; l<d<L t 



2(1 
d 



(2.2.17) 



where d\L t means that d divides L t and <fr(n) is the Euler function, equal to the number of 
positive integers not exceeding the positive integer n and relatively prime to n. 

Substituting eq. ( [2.2. 17|) into eq. (|2.2.14j) then yields a conjecture for N z , sq ,pF,L t (again 
for arbitrary L t ), namely 



N: 



Z,sq,PF,L t 



(it/2) + N Z,tri,PF,L t 



for even L t 
for odd L t . 



(2.2.18) 



where N z ,tri,PF,L t is given by eq. ( p.2.17|) . 



Note that for each strip graphs of a given lattice A with some specified boundary con- 
ditions BC for which one has results from ||18|| ) |l7j and the present work, the ratio of the 



number of A's in the chromatic polynomial divided by the number of A's in the full Potts 
model partition function decreases exponentially fast for L t and vanishes in the limit of large 

hm N J^ BC ^ = . (2.2.19) 

Lt^oo N z ,A,BC,L t 
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3 Potts Model Partition Functions for Square- lattice 
Strips with Free Boundary Conditions 

The Potts model partition function Z(G,q,v), equivalent to the Tutte polynomial, for a 
square-lattice strip of width L t and length L £ = m vertices with free boundary conditions is 
given by 

%,) F x(L t)F (M = v^-H-T^.u (3.1) 
= w T -T m " 1 -M (3.2) 

where w T = v 5 - ■ H . (No confusion should result between the variable v and the vector v.) 
Here T is the transfer matrix, and H is the part of the transfer matrix corresponding to the 
edges (bonds) lying within a single layer. This partition function was calculated for arbitrary 
q, v, and Le with L t = 2 in |T^] and with L t = 3 in |13| . It was also calculated for arbitrary q 
and v for large finite strips in ||19|| . We review the calculations for arbitrary Li with L t = 2, 3 



here in our present notation. We have also calculated the Potts model partition function for 
arbitrary Li for the strip with widths L t = 4 and L t = 5 having free boundary conditions 
and have used these for the analyses of partition function zeros given later in this paper, but 
the results are too lengthy to list here. They are available from the authors on request and 
in the mathematica file transf er_tutte_sq.m that is available with the electronic version 



of this paper in the cond-mat archive at http : / /www . lanl . gov 



3.1 L t 



The number of elements in the basis is trivially equal to two: P = {l,^^}. In this basis, 
the transfer matrices and the other relevant quantities are given by 



T 
H 



w 



u 



q 2 + 3qv + 3v 2 (l + v)(q + 2v) 
v 3 v 2 (l + v) 

1 

V 1 + v 

q 

q(q + v) 
q(l + v) 

1 




(3.1.1) 
(3.1.2) 
(3.1.3) 
(3.1.4) 
(3.1.5) 



In [59, 3(]] a generating function method was presented for chromatic polynomials and 
was generalized to the full Potts model partition function in [0. For a strip graph of type 
G, including specification of width and boundary conditions, with arbitrary length Li = m, 
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and denoting the particular member of this class of graphs with length m as G m , one has 



T(G, q, v,z) = ^T t z ( G m, q, v)z r 



(3.1.6) 



m=0 



The generating function is a rational function of the form 

Af{G,q,v,z) 



with 



and 



T(G,q,v,z) 



V(G,q,v,z) 



3=0 



{q,v)z 3 



(3.1.7) 



(3.1.8) 



V(G,q,v,z) = l + Y^baM^)* (3-1.9) 

3=1 

where the A G j and b G j are polynomials in q and v that depend on the type of strip (but 
not its length) of respective degrees and dp in the auxiliary expansion variable z. The 
factorization of the denominator yields the A's for a given type of strip: 



V(G,q,v,z) 



d-D 
3=1 



\ G)j (q,v)z) 



(3.1.10) 



In [5"9|-||12|| a convention was used that allowed for a more general (inhomogeneous) form of 
strip graph J(Y[ H)I in which the end subgraphs / and J could be different from the subgraph 
H that is repeated to form the strip. Since this generality is not needed for homogeneous 
lattice strip graphs, a convention was introduced in [15j according to which the m = 
term in the expansion (|3.1.6|) is taken to be the line graph Ti t rather than a column of 
squares of height L t . This simplifies the A G j polynomials in Af(G, q,v,z) (the denominator 
T>(G,q,v, z) is independent of this convention). Here we shall use the convention of ||15|| . 
Further, to shorten the notation, we shall suppress the G-dependence on TV, V, the b G ,j, and 
A G j where this is obvious. Then for the present strip graph of width L t = ' 



V 

M 



where JT2] 



A 
A l 



1 + b x z + b 2 z 2 
A + A x z 



(v 3 + 4f 2 + 3qv + q 2 ) 



v 2 (l + v)(q + v) 2 

v(q + v) 
-q 2 v 2 (l + v) 



(3.1.11) 
(3.1.12) 



(3.1.13) 
(3.1.14) 

(3.1.15) 
(3.1.16) 
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In this case we have a four- dimensional basis given by P = {1, Si^ + #2,3, ^1,3, Si. 
basis the transfer matrix is given by 



/ 



T 



11 



T 



12 



T 



13 



v 3 (q + 2v) v 2 {l + v) (q + 3v) 
v 4 2v 3 (l+v) 
v 5 2v 4 (l+v) 



,2/1 1 „,\2 



v 3 (2 + v) v 2 (l + v) 
v 2 (q + 3v + v 2 ) v 2 (l + v) 2 
v 4 (2 + v) v 3 (l + v) 2 J 



where 



Tis 



(q + 2v)(q 2 + 3qv + Av 2 ) 
2(1 + v)(q 2 + Aqv + 5v 2 ) 
q 2 + 5qv + 8v 2 + qv 2 + 3v 3 



(l+v) 2 (q + 3v) 
The rest of the matrices and vectors are given by 



H = 



( 1 

v 






1 + v 




w = 



( 

2q 2 

\ 1 J 
( 

2q( 
Q( 

/1\ 




In the generating-function formalism 







1 

\v 2 2v(l + v) v(2 + v) (l + v) 2 J 

3 \ 



q- 
2q 2 

q 

q(q + v) 2 
2q(l + v)(q + v) 
q(q + 2v + v 2 ) 



J 



u 



V = 1 + b x z + b 2 z 2 + b 3 z 3 + b 4 z 4 
N = A + Aiz + A 2 z 2 + A 3 z 3 



where 



61 = 



5q 2 v - 12qv 2 



15v< 



qv" 



Qv 4 
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b 2 = v 2 {q + v)(2q 3 + 13q 2 v + q 3 v + 30qv 2 + 8q 2 v 2 + 32v 3 + 21qv 3 + q 2 v 3 

+30v 4 + Aqv A + 10v 5 + v 6 ) (3.2.13) 

63 = -v 4 (l + v)(q + v) 2 (q 3 + 9q 2 v + 19qv 2 + 4gV + 15v 3 + 10qv 3 + q 2 v 3 

+10v 4 + 2qv 4 + 2t> 5 ) (3.2.14) 

6 4 = v 7 (l + v) 3 (q + v) 5 (3.2.15) 

A = q(q + v) 2 (3.2.16) 

A x = -qv 2 (-v 4 + Qq 2 v 2 + q 2 v 3 + q 3 v + 2qv 4 + 7qv 3 + 9q 2 v + 9qv 2 + 2q 3 ) (3.2.17) 

A 2 = qv 4 (l + v)(q + v)(2qv 3 -v 3 + q 2 v 3 + 4qv 2 + 4q 2 v 2 + 7q 2 v + q 3 ) (3.2.18) 

A 3 = -q 3 v 7 (l + v) 3 (q + v) 2 (3.2.19) 

4 Potts Model Partition Functions for Strips with Cylin- 
drical Boundary Conditions 

The Potts model partition function for a square-lattice strip of fixed width L t vertices and 
arbitrarily great length L^ = m vertices, with cylindrical boundary conditions, is given by 

Z(L t ) P .(L e )M,v) = v T -H.T m - l .u (4.1) 
= w T ■ T m ~ l ■ u (4.2) 

where w T = ■ H. Here we give results for the cases L t = 2 and L t = 3. We have 
also calculated this partition function for L t = 4, 5, but the results are too lengthy to 
list here. They are available from the authors on request and in the mat hematic a file 
transf er_tutte_sq.m that is available with the electronic version of this paper in the 



cond-mat archive at |http : //www . lanl . gov 



4.1 L t = 2 

The square-lattice strip of width 2 P is equivalent to the square-lattice strip of width 2 F when 
v — — 1 (chromatic polynomial), but this equivalence does not hold for nonzero temperature. 
Furthermore, the square-lattice strip of width 2p with coupling v is easily seen to be equiva- 
lent to an inhomogeneous square-lattice strip of width 2p with coupling v in the longitudinal 
direction and v' — v(2 + v) in the transverse direction. 

The number of elements in the basis is two: P = {1, ^1,2}- In this basis, the matrices and 
vectors are given by 

/ q 2 + Aqv + 5v 2 + qv 2 + 2v 3 (1 + v) 2 {q + 2v) \ . . 

1 ~ \ v 3 (2 + v) v 2 (l + v) 2 J l4 - i " ij 

H = \v{2 + v) (l + v) 2 ) = {v {l + v)J (4L2) 
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q 2 



(4.1.3) 



w = ( ^ + 2v \f) ) (4.1.4) 



Q 

q{l + vf 

a = (J) (4.1.5) 

In the generating-f unction language we have 

V = l + b lZ + b 2 z 2 (4.1.6) 
AT = A + A lZ (4.1.7) 



where 



61 = -{q 2 + Aqv + Qv 2 + qv 2 + Av 3 + v A ) (4.1.8) 

6 2 = v 2 {l + vf{q + v) 2 (4.1.9) 

= g(g + w 2 + 2w) (4.1.10) 

A 1 = - q 2 v \l+v) 2 (4.1.11) 

4.2 L t = 3 

For the strip with width L t = 3 and cylindrical boundary conditions (PF3) there are three 
elements in our basis: P = {l,^i,2 + ^2,3 + ^1,3,^1,2,3}- In this basis the transfer matrix is 
given by 

(Tn T12 T13 

v 3 (q + Av + v 2 ) v 2 {l + v)(q + 7v + 3v 2 ) v 2 {l + vf | (4.2.1) 
v 5 (3 + v) 3v 4 (l + v)(2 + v) v 3 (l + v) 3 

where 

Tn = q 3 + 6q 2 v + 15qv 2 + 16v 3 + qv 3 + 3v 4 (4.2.2) 

T 12 = 3(l + v)(q 2 + 5qv + 8v 2 + qv 2 + 3v 3 ) (4.2.3) 

T 13 = (l+v) 3 (q + 3v) (4.2.4) 

The rest of the matrices and vectors are given by 

1 

H = I w 1 + v 1 (4.2.5) 
v 2 (3 + t;) 3w(l + w)(2 + t;) (1 + v) 3 

g 3 

u = ( 3q 2 j (4.2.6) 
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q{q 2 + 3qv + 3v 2 + v 3 ) 

US = \ 3 q (l + v )( q + 2v + v 2 ) | (4.2.7) 
q(l + v) 3 

n = I o I (1.2.8) 



In the generating-f unction language we have 



V = 1 + hz + b 2 z 2 + b 3 z 3 (4.2.9) 
J\f = A + A lZ + A 2 z 2 (4.2.10) 



where 



b x = -(q 3 + 6q 2 v + I6qv 2 + 2Av 3 + 2qv 3 + 16f 4 + 6t> 5 + v 6 ) (4.2.11) 
b 2 = v 2 (l + v){q + v){q 3 + 10q 2 v + 26qv 2 + 5q 2 v 2 + 2Av 3 + 20qv 3 + q 2 v 3 + 26v 4 

+5qv 4 + 10v 5 + v 6 ) (4.2.12) 

6 3 = -v 5 (l + v) 4 (q + v) 4 (4.2.13) 

A = q(3v 2 + 3qv + v 3 + q 2 ) (4.2.14) 

A x = -qv 2 (l + v)(3qv i - 2v 4 + q 2 v 3 + 10qv 3 - 3v 3 + 5q 2 v 2 + 9qv 2 

+8q 2 v + q 3 ) (4.2.15) 

A 2 = q 3 v 5 (l + v) 4 (q + v) (4.2.16) 



5 Partition Function Zeros in the q Plane 

In this section we shall present results for zeros and the continuous accumulation set (singular 
locus for the free energy) B in the q plane for the partition function of the Potts antiferromag- 
net on square-lattice strips with free longitudinal and free or periodic transverse boundary 
conditions. 



5.1 Free Transverse Boundary Conditions 

The L t = 2 strip with free boundary conditions was previously studied in []T2"|. It will be 



useful to review some of these results. In Fig. |l] we show the accumulation set B formed 
by the partition function zeros in the q plane for the Potts antiferromagnet on the free 
strip of width L t = 2 in the limit of infinite length. For comparison we show the partition 
function zeros calculated for a finite strip with length = 20. Evidently, except for obvious 
discrete zeros such as the zero always present at q — 0, these zeros lie close to the asymptotic 
curves. For v = —1, the continuous locus B = 0, with the zeros accumulating at the discrete 
complex-conjugate points q = e t7T ^ 3 . As v increases above —1, i.e. the temperature for the 
Potts antiferromagnet increases above zero, B forms complex-conjugate arcs whose endpoints 
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nearest to the real axis approach this axis. As was noted in fp72| , when a = e K = v + 1 
increases through the value (3/4) 2 , i.e. v = —7/16, these endpoints touch the real axis, and 
in the interval —7/16 < v < 0, B consists of the above-mentioned complex-conjugate arcs 
together with a real line segment. As v — > CT, the locus B shrinks and finally degenerates 
to the point at q = for v = (see Figs. 3,4 in [1^]). The ferromagnetic case 1 < v < oo 



has been studied in 12 and will not be discussed here. 



in 



Partition functions and zeros in the q plane for finite temperature were previously studied 
|19|1 for the free L t = 3 and L t = 4 strips. In Fig. || we show B in the q plane for the 
infinite-length limit of the L t = 3 free strip and, for comparison, partition function zeros 
calculated for a finite strip with length Lg = 30, for the Potts antiferromagnet. The zeros in 
this figure are similar to those in Fig. 3.3 in [jnj. For zero temperature, i.e., v = —1, B is 
comprised of the union of a self-conjugate arc passing through q = 2 and a complex-conjugate 
pair of arcs (see Fig. 3(a) in J5U|). As the temperature increases above zero, the arcs tend 
to move together and contract toward the origin. In Fig. |] we show B for the infinite-length 
limit of the free L t = 4 strip and, for comparison, partition function zeros calculated for 
a finite strip with length Li = 40, for the Potts antiferromagnet. The zeros in this figure 



are similar to those in Fig. 3.9 in |lSf| . For zero temperature, B consists of several arcs 



together with a small real line segment (see Fig. 3(b) in |59f). There are complex-conjugate 



triple points evident on this locus. Such triple points were studied in in the context of 
complex-temperature phase diagrams, and it was shown how they arise when curves on the 
singular locus B meet in such a manner that, as one travels along given curve(s) past the 
intersection point the eigenvalues whose degeneracy in magnitude defines the curve cease 
to be dominant eigenvalues (see Figs. 1 and 2 of fl4"0"f ). The nature of these triple points 
was further analyzed in (T7]]. If one considers B as the union of the various curves and line 
segments that comprise it, then a triple point is a multiple point (intersection point) on B, 
since it lies on multiple branches of B. In a different nomenclature, if one considers each 
of the algebraic curves comprising B individually, including the portions where the pairs 
of degenerate-magnitude A's are not dominant so that these portions are not on B, then 
such triple points are not multiple points on each individual algebraic curve, since these 



individual curves pass through the triple point as shown in Fig. 2 in p0[ . As was found for 
the strips with L t = 2, 3, as the temperature increases, the curves forming B contract toward 
the origin. One can observe, in particular, that for sufficiently high temperature the triple 
points disappear. 

We proceed to the free strip of width L t = 5. In Fig. |] we show B for the infinite-length 
limit of the free L t = 5 strip and, for comparison, partition function zeros calculated for a 
finite strip with length Lg = 40, for the Potts antiferromagnet. The special case v = — 1 
(chromatic polynomial) was studied in |]72| , [17| . The results exhibit the general trend that 
for v = — 1 as Lt increases, the arcs on B move together, and the endpoints nearest to 
the origin approach this point. Just as was seen for narrower widths L t , as the temperature 
increases, the arcs tend to come together further, the prongs protruding to the right contract 
and eventually disappear, and the entire locus contracts toward the origin. While B has an 
oval-like shape for v — — 1, it becomes more round as v — > 0~. 
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5.2 Periodic Transverse Boundary Conditions 



We next consider the strips with cylindrical boundary conditions. In Figs. [5|-[8] we show B in 
the q plane for the infinite-length limits of the cylindrical strips with widths L t = 2 through 
L t = 5 and partition function zeros on long finite strips of the respective widths, for the 
Potts antiferromagnet. For the v — — 1 special case with L t = 2, 3, B degenerates to discrete 
points. The v — — 1 special case for L t = 4 was previously studied in []60| and for L t = 5 in 



|72] , |TT|] (as well as higher widths in these papers, going up to L t = 7 in |I7| and L t = 13 in 
75fl). The cases L t = 3 and L t = 4 for arbitrary temperature were studied in ]l^] and our 
zeros for L t = 4 are similar to those in Fig. 3.12 of Our present results show that the 
locus B tends to become somewhat more complicated as L t increases. 



6 Partition Function Zeros in the v Plane 

We next present results for complex-temperature (Fisher) zeros and the continuous accu- 
mulation set (singular locus for the free energy) B in the v plane for the partition function 
of the Potts model on square-lattice strips with free longitudinal boundary conditions and 
free or periodic transverse boundary conditions. Since the infinite-length limits of the strips 
considered here are effectively one-dimensional systems and since a one-dimensional spin 
system with short-range interactions does not have any finite-temperature phase transition, 
it follows that the locus B where the free energy is singular does not intersect the real f-axis 
in the interval — 1 < v < oo corresponding to nonzero temperature. Hence there cannot be 
any physical finite-temperature phase with ferromagnetic or antiferromagnet ordering, and 
the complex-temperature phase diagram involves only the (complex-temperature extension 



of the) paramagnetic phase, together with possible O phases, in the nomenclature of In 
Fig. H we show the Fisher zeros for long finite strips with free boundary conditions, and the 
approximate respective loci B for the infinite-length limits of these strips with L t ranging 
from 2 through 5. (In Fig. ^](d) the locus B is shown only for the physical range Re(t> ) > — 1.) 
Corresponding results are presented in Fig. |10| for strips with cylindrical boundary condi- 
tions. For each specific width and type of boundary condition we show curves and zeros 
for q — 2, 3, and 4. Although some of these plots are rather complicated, the reader can 
distinguish the curves for different values of q by their close association with the zeros, which 
are labelled with different symbols, namely squares, circles, and triangles for q = 2, q = 3, 
and 9 = 4, respectively. Furthermore, in the postscript files for the figures in the cond-mat 
archive, the curves for different q values have different colors: black for q = 2, red for q = 3, 
and green for q = 4, thereby rendering them easily distinguishable. To complement the types 
of comparisons that can be made with these plots, we show in Fig. [11] and the loci B and 
zeros for various widths plotted together for each value of q from 2 to 4. 

There are a number of interesting features that are evident in these plots. In cases where 
one uses boundary conditions that are self-dual, a subset of the Fisher zeros lie exactly on 



the circle \v\ = ^fq [^3|],|24|, ^0), [77| . It has also been proved that in the limit q — > oo 
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B consists of the unit circle in the complex ( plane, where 



c = — 



(6.1; 



Furthermore, the approach to this limit has been studied using exact solutions for Z(G, q, v) 
on infinite-length, finite- width self-dual strips in [[77]]. Although the open and cylindrical 
strip graphs with finite Li used here are not self-dual, in the case of cylindrical boundary 
conditions, the effect of this non-self-duality is significantly reduced in the limit Li — > 00. 
Indeed, we find that a subset of the Fisher zeros lie on the respective circles \v\ = ^Jq, as 
is evident in the figures. In [23. 24 , 26, 19. 2C] these zeros were studied for finite slices of 



the two-dimensional square lattice, and the resultant Fisher zeros form a circular pattern in 
the v plane which one may infer would become, in the thermodynamic limit, the right-hand 
boundary of the (complex-temperature extension of the) paramagnetic phase, separating it 
on the right from the analogous extension of the ferromagnetic phase. Since the Potts model 
does not have a ferromagnetic phase on the infinite-length finite-width self-dual strip graphs 
used in J77J , what was found there was that the arcs on B in the v plane end at certain angles 
that depend on L t and q. Let us introduce polar coordinates for the variable £ defined in 
(|6.1D: ( = \C\e ld and denote the value of 9 at the right-hand endpoint of the arc of B with 
\(\ = 1 and lm(v) > as 9 ae . For fixed L t , the endpoint angle 9 ae decreases with increasing 
q, and for fixed q, 9 ae decreases with increasing L t . Further, for fixed L t , lim 9 _»oo0 ae = 0, 
and for fixed q, lim^-^ 9 ae = p5[ [77]. For our present strips, we find the following results 
showing the same quantitative trends. For each strip we list three approximate arc endpoint 
angles corresponding to q = 2, q = 3, and q = 4: for 2 P x oo F , 9 ae ~ 74°, 65°, and 60°; for 
3 P x oo F , 9 ae ~ 50°, 40°, and 35°; for 4 P x oo F , 9 ae ~ 38°, 29°, and 24°; and for 5 P x oo F , 
9 ae ~ 31°, 23°, and 18°. 

For q = 2 in the figure comparing different widths for cylindrical boundary conditions, 
Fig. 0(a), one can see that the zeros lie not just on the circle \v\ = v2 but also on the 
other of the Fisher circles, \v + 2| = ^/2, and one can observe how the endpoints on these 
circles move closer to the real axis as L t increases. As L t — > 00, these would then close to 
produce the full Fisher circles for the model on the infinite square lattice. Finally, for the 
values q = 3,4, one can observe prongs on the left-most part of B that are reminiscent of 
analogous prongs whose endpoints were determined from analyses of low-temperature series 



expansions for the infinite 2D square lattice in |f26|. 



7 Internal Energy and Specific Heat 

From the partition function and the resultant reduced free energy for the finite and infinite- 
length strip, ( |1.4j ) and ( |1.5|) , it is straightforward to compute the (internal) energy E and 
specific heat C as 
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and 



C 



dE 
df 



k B K 2 (v + l) 



df | | d 2 f 
dv dv 2 



(7.2) 



In the limit Lg — > oo, since only the dominant eigenvalue X d of the transfer matrix contributes 
to the free energy, one has 



E 



f = -r^x d 

J(v + 1) d\ d 
L t X d dv 



and 



C 



k B K 2 {v + l) 
LtXd 



(v + 1) 



d 2 X d (v + 1) 
dv 2 



X n 



dXd 
dv 



+ 



dXd 
dv 



(7.3) 
(7.4) 

(7.5) 



As noted, since the infinite-length limits of the strips considered here are quasi-one- dimensional 
systems with free energies that are analytic for all finite temperatures, it follows that the 
dominant eigenvalue is the same on the whole semi-axis Im(u) = 0, Re(t> ) > —1. Further- 
more, as discussed in ||12| , the free energy and thermodynamic functions such as the internal 
energy and specific heat for infinite-length finite width lattice strips are independent of the 
longitudinal boundary conditions (although they depend on the transverse boundary condi- 
tions). 

For convenience we define a dimensionless internal energy 



E r 



E 
1 



(7.6) 



Note that sgn(E r ) is (i) opposite to sgn(E) in the ferromagnetic case where J > for which 
the physical region is < v < oo and (ii) the same as sgn(E) in the antiferromagnet case J < 
for which the physical region is — 1 < v < 0. We recall the high-temperature (equivalently, 
small-|i^|) expansion for an infinite lattice of dimensionality d > 2 with coordination number 

~- + ^^ + 0(/d . (7.7) 



E 
1 



A 

2 



q q* 

In passing, it should be noted that in papers on the q = 2 Ising special case, the Hamiltonian 
is usually defined as Hi = —Ji Yliij) a i a j w hh <jj = ±1 rather than the Potts model definition 
( |1 . 2|) ; the isomorphism between these conventions involves the rescaling 2Kj = K, where 
Kj = (3Ji. Furthermore, Ej = — Jfaaj) rather than the Potts definition E = —J{5 aiaj ), 
where (ij) are adjacent vertices. Hence, for example, for q = 2, with the usual Ising model 
definitions, Ej(v = 0) = rather than E = — JA/(2g) and the high-temperature expansion 
is £ 7 = -J(A/2)[K + 0(K 3 )] rather than the q = 2 form of (Q. 
We also define the reduced function 



C 



H 



c 

knK 2 



(7- 
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We show in Figures [T3Hl5| the reduced energy E r and the function Ch that enters in the 
specific heat for q = 2 through q = 4 for the square-lattice strips of widths 2 < L t < 5 with 
free and cylindrical boundary conditions. We recall that the Potts model on the infinite 
square lattice has a phase transition from the paramagnetic high-temperature phase in the 
interval < v < v p ^fm to the low-temperature phase with ferromagnetic long-range order in 
the interval v Pt pu < v < oo, where 

Vp,FM = y/q ■ (7.9) 

This phase transition is continous (second-order) for q < 4 and first-order for q > 4. In the 
figures plotting the function Ch entering in the specific heat one can see how the maxima at 
the respective values of v ~ v p ^m for various q increase as the width of the strip L t increases. 
In the limit L t — > oo, these maxima diverge, corresponding to the divergence of the specific 
heat at v p ^fm- Since the free strips and the cylindrical strips with even L t are bipartite graphs, 
there is a well-known isomorphism between the ferromagnetic and antiferromagnetic Ising 
models, and one sees corresponding maxima at the values of v obtained by the replacement 
K -> -K. 

A general property for E r calculated for the infinite-length limit of the lattice strips with 
free transverse boundary conditions is that E = for v = —1, i.e., the zero-temperature 
limit of the Potts antiferromagnet. This is evident in the figures. For (the infinite- length 
limit of) the square lattice strips with periodic transverse boundary conditions, the value of 
E = at v = — 1 for q > x, the chromatic number for these strips, which is 2 if L t is even 
and 3 if L t is odd. For q = 2 and L t odd, there is frustration, and hence E is larger than 
zero at v = —1. 

In the limit of infinite temperature on the infinite square lattice, E r = 2/q, as is evident 
from the expansion ( |7.7| ). The infinite- length strips with cylindrical boundary conditions 
have uniform coordination number A = 4, and one sees the agreement with the formula 
E r = 2/q in the figures. The infinite- length strips with open boundary conditions do not 
have a uniform coordination number (this is equal to 3 for the vertices on the upper and 
lower sides and 4 for the vertices in the interior of the strip). However, one can see the 
approach to the above formula as L t increases. 

For the g-state Potts ferromagnet on the infinite-length finite- width strips with cylindrical 
boundary conditions, we find that the curves for the energy E cross at a unique value of 
v depending on q but independent of the width L t ; furthermore, this value is equal to the 
value v = v Pt FM at which the phase transition occurs on the infinite square lattice. Let us 
denote 

E p = E(v = v p , FM ) (7.10) 

for the internal energy of the infinite-length finite-width strips evaluated at the value of v 
given in ( fT.9|) . A careful evaluation of the value of E p reveals that for the range of q < 4 where 
the Potts ferromagnet on the square lattice has a continuous second-order phase transition, 
this value is independent of the width L t and is equal to the value for the infinite square 
lattice at the ferromagnetic critical point, given by |^7| 

£7 C = -J(1 + <T 1/2 ) • (7.11) 
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This behavior could have been anticipated in the Ising case, q = 2. Finite-size relations for 
statistical mechanical models using methods of conformal field theory have been discussed in 
88|-f92]. For the Ising model, the asymptotic expansion of the internal energy evaluated at 



v = v Pi fm on a torus of length Le and width L t takes the following form (where we shall use 
the subscript c to denote the fact that this would be the critical value in two-dimensional 



thermodynamic limit) [P3l , fl88| , |94j 195] 



£ c (W) = £ c + £§±^> (7.12) 

where 



i=0 



P = T (7-13) 
is the aspect ratio of the torus. The infinite- length torus corresponds to the limit p —>■ oo. 



It has been shown [54, [5^] that in this limit 

£i(oo) = E 3 (oo) = £5(00) = . (7.14) 

Furthermore, one can see that all the correction terms have a factor proportional to 626364 
where 9i is the usual ^-function 6i(z,r) (e.g., [Q) evaluated at z = 0, where r = ip, i.e. 
q nome = e tT = e _7rp (see, e.g., Appendix B of p5|). In the limit p — > 00, from the basic 
definitions of these theta functions, it follows that 62 — > and 63,64 — > 1. Hence, for the 
Ising model on an infinitely long cylinder 

E c (L t ,oo)=E c (7.15) 

for all widths L t . This (together with the property that the thermodynamic functions are 
independent of the longitudinal boundary conditions in the limit of infinite length) explains 
why all of the curves for E r on the infinite-length, finite-width strips with cylindrical bound- 
ary conditions cross at the point v = v Pt pM where the the energy attains the critical value 
on the infinite square lattice, E c . Our numerical results suggest the inference that all the 
finite-size corrections to the internal energy at the respective ferromagnetic critical points 
also vanish (for an infinitely long cylinder) for the other two values of q, i.e., q = 3 and 
q = 4, where the Potts ferromagnet has a second-order phase transition on the infinite 
square lattice. We have also checked and confirmed that the value of the internal energy 
at this crossing, (E r ) c (v = v p ^m) is equal to 1 + g -1 / 2 for q > 4, where the model on the 
infinite square lattice has a first-order transition, so that the limits of the internal energy at 
v = Vp.FM are different when approached from high and low temperature P7[ . 

Given the bipartite property of the strip graphs with cylindrical boundary conditions 
and even L t and the consequent equivalence of the Ising ferromagnet and antiferro magnet, 
and taking into account our discussion above, it follows that for q = 2 the internal energy 
curves should exhibit a crossing in the antiferromagnetic regime for even L t , and our results 
agree with this. The point where this occurs is v = v Pi afm = V% — 2, and at this point the 
internal energy is precisely equal to the value on the infinite square lattice, E r = l- l/y/2. 
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8 Appendix 

In this appendix we recall the connection between the Potts model partition function Z(G, q, v) 
and the Tutte (also called Tutte/Whitney) polynomial T(G,x,y) for a graph G = G(V,E), 
given by §-0 

T(G,x,y)=J2(x-l) k{G ' ) - k{G) (y-V C{G ' ) (8-1) 

G'CG 

where G' is a spanning subgraph of G, and k(G'), e(G'), and n(G') = n(G) denote the 
number of components, edges, and vertices of G', and 

c(G') = e{G') + k{G') - n{G') (8.2) 

is the number of independent circuits in G' . As stated in the text, k(G) = 1 for the graphs 
of interest here. Now let 

x=l + - (8.3) 

v 

and 

y = a = v + 1 (8.4) 

so that 

g = (z-l)(y-l) . (8.5) 

Then 

Z(G,q,v) = (x-l)^ G \y-l) n ^T(G,x,y) . (8.6) 
For a planar graph G the Tutte polynomial satisfies the duality relation 

T(G,x,y) = T(G*,y,x) (8.7) 

where G* is the (planar) dual to G. 

There are several special cases of the Tutte polynomial that are of interest. One that 
we have commented on in the text and in previous papers is the chromatic polynomial 
P(G,q). This is obtained by setting y = 0, i.e., v = —1, so that x = 1 — q; the cor- 
respondence is P(G,q) = (—q) k ^ G \—l) n ^T(G,l — q,0). A second special case is the 
flow polynomial || |^, |TD| F(G,q), obtained by setting x = and y = 1 — q: F(G,q) = 
(_iy{G)-n(G)+KCf)T(G, 0, 1 - q) For planar G, given the relation (|8"7?|) , the flow polyno- 
mial is, up to a power of q, proportional to the chromatic polynomial on the dual graph 
G: F(G,q) oc P(G*,q). A third special case for x = 1 is the reliability polynomial (e.g., 
||). Consider a connected graph G(V,E) and now let each edge of G be present with 
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probability p (and hence absent with probability 1 — p); then the probability that exists a 
path connecting any two vertices in G is given by the (all-terminal) reliability polynomial 
R(G,p) = J2iicGP e ^(^ ~ p) e ( G >~ e ( H \ where H denotes a connected spanning subgraph of 
G. In turn, this is expressed as a special case of the Tutte polynomial according to 

R(G,p) =p n(G) - 1 (l-j 5 ) e(G) - n(G)+1 T(G,l,(l-p)- 1 ) (8.8) 

Thus, our results for Z(G,q,v) and hence, via (|8.6|) , for T(G,x,y) for these lattice strips 
yield, as a special case, the reliability polynomials for the strips. A conjecture for the complex 
zeros of the reliability polynomial was given in |)7| (see |67|] for a recent discussion). 



For a given graph G = (V,E), at certain special values of the arguments x and y, the 



Tutte polynomial T(G,x,y) yields quantities of basic graph-theoretic interest 0-JlO[- We 
recall some definitions: a spanning subgraph was defined at the beginning of the paper; a 
tree is a connected graph with no cycles; a forest is a graph containing one or more trees; 
and a spanning tree is a spanning subgraph that is a tree. We recall that the graphs G that 
we consider are connected. Then the number of spanning trees of G, Nst{G), is 

N ST (G) = T(G, 1,1) , (8.9) 

the number of spanning forests of G, N S f(G), is 

N SF (G)=T(G,2,1) , (8.10) 

the number of connected spanning subgraphs of G, N GSSG (G), is 

N GSSG (G) = T(G, 1,2), (8.11) 

and the number of spanning subgraphs of G, Nssg(G), is 

N SSG (G)=T(G, 2,2) . (8.12) 

From the duality relation (|8.7|) , one has, for planar graphs G and their planar duals G* , 

N ST (G) = N ST (G*) , N SSG (G) = N SSG (G*) (8.13) 

and 

N SF (G) = N CSSG (G*) . (8.14) 

In previous works |l2l-||16|| we have given resultant formulas for these graphical quantities 
for the families of graphs considered therein. It is straightforward to do the same for the 



square-lattice strips considered here. For spanning trees on lattice sections, see also [98, |99 . 
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Figure 1: Partition-function zeros in the complex q plane for the 2p x 20f square-lattice 
strip and resultant accumulation set B for the 2f x oop strip (singular locus for the free 
energy) for several values of the temperature-like parameter v: —1 (□), —0.75 (o), —0.50 
(A), —0.25 (O), and —0.10 (©), where the symbols correspond to the zeros calculated for 



the finite strip. The zero at q 



is present for all v. 
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Figure 2: Partition-function zeros in the complex q plane for the 3f x 30f square-lattice 
strip for several values of the temperature-like parameter v: —1 (□), —0.75 (o), —0.50 (A), 
—0.25 (O), and —0.10 (©). The loci B for the 3f x oop strip for these values of v are also 
shown. 
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Zeros sq lattice L. = 4 
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Figure 3: Partition-function zeros in the complex q plane for the for the 4-p x 40p strip for 
several values of the temperature-like parameter v: — 1 (□), —0.75 (o), —0.50 (A), —0.25 
(O), and —0.10 (©). The loci B for the 4p x oo F strip for these values of v are also shown. 
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Zeros sq lattice L t = 5 F 
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Figure 4: Partition-function zeros in the complex q plane for the 5f x 40f square-lattice 
strip for several values of the temperature-like parameter v: —1 (□), —0.75 (o), —0.50 (A), 
—0.25 (O), and —0.10 (©). The loci B for the 5f x oof strip for these values of v are also 
shown. 
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Figure 5: Partition-function zeros in the complex q plane for the 2p x 20f square- lattice 
strip for several values of the temperature-like parameter v: —1 (□), —0.75 (o), —0.50 (A), 
—0.25 (O), and —0.10 (©). The loci B for the 2p x oop strip for these values of v are also 
shown. 
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Figure 6: Partition-function zeros in the complex q plane for the 3p x 30f square-lattice 
strip for several values of the temperature-like parameter v: —1 (□), —0.75 (o), —0.50 (A), 
—0.25 (O), and —0.10 (©). The loci B for the 3p x oo F strip for these values of v are also 
shown. 
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Zeros sq lattice L. = 4 
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Figure 7: Partition-function zeros in the complex q plane for the 4p x 40f square-lattice 
strip for several values of the temperature-like parameter v: —1 (□), —0.75 (o), —0.50 (A), 
—0.25 (O), and —0.10 (©). The loci B for the 4p x oop strip for these values of v are also 
shown. 
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Figure 8: Partition-function zeros in the complex q plane for the 5p x 40f square-lattice 
strip for several values of the temperature-like parameter v: —1 (□), —0.75 (o), —0.50 (A), 
—0.25 (O), and —0.10 (©). The loci B for the 5p x oo F strip for these values of v are also 
shown. 
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Figure 9: Partition function zeros in the complex v plane for several square-lattice strips: 
(a) 2 F x 20p, (b) 3f x 30p, (c) 4 F x 40p, and (d) 5f x 50p- In each plot we show the zeros 
for several values of the parameter q: 2 (□, black), 3 (o, red), and 4 (A, green) and the 
corresponding loci B for these values of q. 
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Figure 10: Partition function zeros in the complex v plane for several square-lattice strips: 
(a) 2 P x 20p, (b) 3p x 30p, (c) 4 P x 40p, and (d) 5p x 50p- In each plot we show the zeros 
for several values of the parameter q: 2 (□, black), 3 (o, red), and 4 (A, green) and the 
corresponding loci B for these values of q. 
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Figure 11: Limiting curves forming the singular locus £>, in the v plane, for the free energy 
of the q = 2 (a) q = 3 (b), and q = 4 (c) Potts model on the (L t )F x oof square-lattice strips 
with the following widths L t : 2 (black), 3 (red), 4 (green), and 5 (blue). We also show the 
partition-function zeros for the strips {L t )p x (10L t )p for the same values of L t : 2 (□, black), 
3 (o, red), 4 (A, green), and 5 (o, blue). 
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Figure 12: Limiting curves forming the singular locus £>, in the v plane, for the free energy 
of the q = 2 (a) q = 3 (b), and q = 4 (c) Potts model on the (L t )p x oop square-lattice strips 
with the following widths L t : 2 (black), 3 (red), 4 (green), and 5 (blue). We also show the 
partition-function zeros for the strips (L t )p x (10L t )p for the same values of L t . The color 



code is as in Figure [11 
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Energy sq lattice F bc's q=2 Energy sq lattice P bc's q=2 




Figure 13: Reduced internal energy E r = —E/J and specific heat C# as functions of 
the temperature-like parameter v for the q = 2 Potts model on square-lattice strips of size 
(L t )p x oo F an d {Lt)p x oof with 2 < L t < 5. The plot includes both the ferromagnetic 
and antiferromagnetic Potts models, for which the temperature ranges are < v < oo and 
— 1 < v < 0, respectively. 
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Figure 14: Reduced internal energy E r = —E/J and specific heat C# as functions of v 
for the q = 3 Potts model on square-lattice strips of size (L t ) F x oop and (L t )p x oof with 
2 < L t < 5. Notation is as in Fig. O. 
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Figure 15: Reduced internal energy E r = —E/J and specific heat C# as functions of v 
for the q = 4 Potts model on square-lattice strips of size (L t ) F x oof and (L t )p x oop with 
2 < L t < 5. Notation is as in Fig. O. 
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